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Responses to chemical potential of the pseudoscalar meson screening mass and the chiral condensate in lattice 
QCD are investigated. On a 16 x 8 2 x4 lattice with two flavors of staggered quarks the first and second responses 
below and above T c are evaluated. Different behavior in the low and the high temperature phases are observed, 
which may be explained as a consequence of the chiral symmetry breaking and restoration. 
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1. Introduction 

The study of finite baryon density at finite tem- 
perature is not only very useful to understand 
the phase transition between hadron and quark 
gluon plasma, but also quite important for heavy 
ion collision experiments which require theoret- 
ical understanding of hadronic properties at fi- 
nite baryon density and temperature [Q. As the 
fermionic determinant at finite chemical potential 
is complex on the lattice, numerical simulations 
are very difficult but the quenched approximation 
at finite chemical potential can give physically 
misleading results simulations with dynami- 
cal fermions are therefore essential to extract the 
relevant physics j3|. 

We propose a new technique to investigate non- 
zero chemical potential using lattice QCD simu- 
lations. A Taylor expansion in is used in the 
vicinity of zero /i at finite temperature for eval- 
uating masses and chiral condensate. Here we 
obtain the Taylor coefficients directly in a simu- 
lation at fj, = Q, by measuring the derivatives of 
the relevant observables. There is in fact much 



interesting physical information which can be ex- 
tracted from the behavior of a system at small 
chemical potential. Our preliminary results have 
been reported in [Q. 

2. Lattice formulation 

Assume that the spatial hadron correlator C(x) 
is dominated by a single pole contribution, 
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where M = aM and x 
size in the x-direction. 

We take the first and second derivatives of the 
hadron correlator with respect to p, = afJ, = 
(i/(N t T), where /i is the chemical potential. 
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In this work, we consider the flavor non-singlet 
mesons in QCD with two flavors. The hadron 
correlator is then given by 

(#Wff(0)t) = (G) (4) 

= (Tr [P(Pu)nOrP(Pd)on^] } . 

Here P{p) = D[U;p]~ l {D[U;p} is the Dirac op- 
erator) is the quark propagator at finite chemical 
potential, and T is the Dirac matrix which speci- 
fies the spin of the meson. In this paper we study 
the response to the isoscalar and isovector chem- 
ical potentials /is , py by setting 



Ms — Pu — ; Pv — pu 



(5) 



3. Numerical Simulations and Results 

The simulations have been performed at finite 
temperature T/T c G~ [0.9,1.1] on a 16 x 8 2 x 
4 lattice with standard Wilson gauge action and 
with two dynamical flavors of staggered quarks. 
We use the R-algorithm, with quark masses ma = 
0.0125, 0.017 and 0.025. We also use a corner- 
type wall source after Coulomb gauge fixing in 
each (y, z, t)-hyperplanc. 

The first derivative of the pseudoscalar meson 
correlator with respect to the isoscalar chemical 
potential is identically zero. For the isovector 
chemical potential, our simulation values for the 
first derivative are very small in both phases. 

3.1. Response of the pseudoscalar meson 
to the isoscalar chemical potential 

In the low temperature phase, the dependence 
of the mass on ps is small. This behavior is to 
be expected, since, below the critical temperature 
and in the vicinity of zero ps, the pseudoscalar 
meson is still a Goldstone boson. In fact, the 
chiral extrapolation of the isoscalar response is 
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Figure 1. d 2 M 2 / dp 2 s for the pseudoscalar meson 
versus ma at T < T c (J3 = 5.26, triangles) and 
T > T c ((3 — 5.34, circles). Extrapolation to 
ma = is also shown. 



consistent with zero, as shown in Figure This 
is in contrast with the behavior above T c , where 
d 2 M 2 1 dp 2 seems to remain finite even in the chi- 
ral limit. In addition, our results suggest that the 
response of the coupling A (eq.(|l])) is small below 
T c . 

Above T c , we first note that the correlator and 
its response are still well fitted by the single pole 
formulae, Eqs. (|l]-[|). The screening masses are 
manifestly larger than those below T c . As pointed 
out above, the response of the mass above T c be- 
comes large and positive, reflecting the fact that 
the pion is no longer a Goldstone boson and in- 
dicating chiral symmetry restoration. 

3.2. Response of the pseudoscalar meson 
to the isovector chemical potential 

In the presence of the isovector chemical poten- 
tial, 7r + and 7T~ may have different masses. Here 
we consider the ir + ( ud ) meson as in Eq. (|]) . 

An interesting point in this respect is that 
the second derivative of the mass is negative in 
the low temperature phase, in contrast with the 
isoscalar potential case: the mass tends to de- 
crease under the influence of the isovector chem- 
ical potential. This may be explained by the ob- 
servation that, for low temperature and chemical 
potential above the pion mass, a Goldstone mode 
can appear §. The behavior of M is more clearly 
shown by a Taylor expansion: At j3 — 5.26 and 
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Figure 2. Second responses d 2 M/d{ig and 
d 2 M/djly of the pseudoscalar meson mass at 
ma = 0.025. 



ma = 0.017, the data suggest 
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In the high temperature phase, the depen- 
dence of the masses on /iy decreases. Since the 
pseudoscalar meson becomes heavier, the phase 
boundary to the pion condensate phase is further 
away from the [iy = axis. The weaker responses 
may be understood from this point of view. 
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Figure 3. Second responses d 2 {-tptp) / dfi s and 
d 2 (tpip) I dfly at ma = 0.025. Note that the ab- 
scissa of coincident points have been splitted. 



3.3. Responses of the chiral condensate 

We also measure responses of the chiral con- 
densate to the isoscalar and isovector chemical 
potentials. The first responses to both potentials 
are identically zero. Figure [| shows our prelim- 
inary results for the second responses. We have 
found they are small and negative in both phases. 
Near T c there is a large change. The critical tem- 
perature tends to decrease under the influence 
of nsy- Thus, in the low temperature phase, 
turning on the chemical potential brings the sys- 
tem closer to the phase transition where chiral 
symmetry is restored, and decreases the chiral 
condensate. At high temperature, because chi- 
ral symmetry is restored, responses of the chiral 
condensate to the isoscalar and isovector chemical 
potential are small. 

4. Conclusions 

We have developed a framework to study the 
response of hadrons to the chemical potential. 
The dependence of the pseudoscalar meson mass 
on us in the chiral limit is consistent with zero at 
low temperature, reflecting the fact that at small 
/is the pion is still a Goldstone boson. For the 
isovector chemical potential, the ud pseudoscalar 
meson mass tends to decrease as a function of y,y 
at a much stronger rate in the low temperature 
phase. This is consistent with a restoration of chi- 
ral symmetry as the isovector chemical potential 
increases. 
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